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Abstract 



J^ ' We point out that the axino predicted as the supersymmetric partner of the 

0^ ■ axion is a good candidate for the recently proposed sterile neutrino cool dark 

FjJ . matter. The axino mass falls into the right range (^ keV) in the context of 

f^ i gauge mediated supersymmetry breaking. A sizable mixing of the axino with 

active neutrinos arises when R-parity violation is allowed and the resulting 

neutrino masses and mixing accommodate the atmospheric and solar neutrino 

p ■ oscillations simultaneously. 
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In a recent paper |l| , a new cosmological model has been proposed to explain the struc- 
ture formation of the Universe. In this model, dark matter consists of non-thermal sterile 
neutrinos with masses about 100 eV to 10 keV which are produced around the big bang 
nucleosynthesis (BBN) epoch through a resonant active-sterile neutrino conversion in the 
presence of a net lepton number asymmetry. Such non-thermal neutrinos are called "cool" 
since their free streaming scale is larger than that of the usual cold dark matter but smaller 
than that of hot or warm dark matter. This cool dark matter (CIDM) model would be a 
good alternative to the conventional cold plus hot dark matter (CHDM) model which has 
been known to be the best cosmological model for the structure formation [Q. In view of 
particle physics theory, the validity of the CHDM model depends crucially on the origin of 
active neutrino masses and mixing. In the CHDM model, 20 % of dark matter is composed 
of hot components which are unarguably taken as three species of active neutrinos. The 
cosmological requirement, fij, ~ 0.2, constrains the sum of three active neutrino masses 
through the relation, J2'^u = 92Qiyh'^ eV, which implies J2''t^u — 4.6 eV, taking h = 0.5. If 
neutrinos take natural hierarchical masses as quarks and charged leptons, the recent atmo- 
spheric neutrino data from the Super-Kamiokande |^ imply Am^^^^ ~ ml^ < 10~^ eV^ and 
thus Y^m,y < 0.1 eV, in which case the CHDM model cannot work. 

The purpose of this paper is to present a particle physics model in which the CIDM 
model is realized in a natural way while neutrinos obtain hierarchical masses which can 
accommodate the atmospheric and solar neutrino data |Q simultaneously. One of popular 
ways to generate neutrino masses is to consider the minimal supersymmetric extension of 
standard model without imposing R-parity in which R-parity and lepton number violating 
bilinear and trilinear couplings give rise to typically hierarchical neutrino masses PJ§]. If 
one introduces in this context the Peccei-Quinn (PQ) mechanism which would be the most 
attractive solution to the strong CP problem ^ , there exits a singlet fermion, called the axino 
a, which is a supersymmetric partner of the Goldstone boson, the axion, of a spontaneously 
broken PQ symmetry. Various cosmological and astrophysical observations are known to 
restrict the scale of PQ symmetry breaking: fa ~ (10^° — 10^^) GeV ||^. As was observed 
first in Ref. @, the axino can be a sterile neutrino, that is, it mixes with active neutrinos 
once R-parity is not imposed. In the context of supergravity [^ where supersymmetry 
breaking is mediated at the Planck scale, the axino mass is quite model-dependent and 



could be in the range, 100 GeV — 1 keV |Tl|. However, in models with gauge mediated 



supersymmetry breaking (GMSB) [T^ where the supersymmetry breaking mediation scale 
is much below the PQ scale fa, the axino mass is determined by the Goldstone nature of the 
axion supermultiplet, and is predicted to be typically in the sub-keV regime in the minimal 
gauge mediation models [T^ . 



We will show that the axino can be a good candidate of cool dark matter as its mass and 
mixing with active neutrinos fall into the right ranges for a reasonable range of parameter 
space in the context under consideration. In particular, the required mixing of the axino 
with active neutrinos can be obtained when one allows the R-parity violating terms whose 
sizes are fixed to explain the atmospheric and solar neutrino masses and mixing. A weak 
point of CIDM scenario would be the requirement of a large lepton asymmetry. It will be 
argued that the demanded lepton asymmetry can be generated through a late-time entropy 
production followed by the R-parity and lepton number violating decays of the ordinary 
lightest supersymmetric particle (LSP), or through the Affleck-Dine mechanism. 



There are various sources for the cosniological axino population. The conventional source 
would be the thermal production. If the reheat temperature Tr after inflation is larger than 
the axino decoupling temperature T^ ~ /„, thermally produced axinos may overdose the 



universe unless the axino mass is less than about 2 keV [0]. Therefore, the axino with mass 
~ keV can be warm dark matter. When Tr < Td, the primordial axinos are inflated away, 
but can be regenerated from thermal background. In this case, the axino warm dark matter 
needs the mass | |15|| ; 



[W^GeVj \ Tr J ^ ^ 

The other possibility is the secondary production of the axinos through the decay of ordinary 



superparticles. In a recent paper |]I6[, it was argued that the axino with mass ~ 10 GeV can 
be cold dark matter when the axinos are produced from the decay of the ordinary neutralino 
LSP. The last one would be the nonthermal production resulting in the axino cool dark 
matter, which we would like to discuss here. 

Let us recapitulate the basic ingredients of the CIDM model realizing the nonthermal 
production ||l|]. The production of a sterile neutrino z/^ (= a in our scheme) which is much 
heavier than active neutrinos can be made by a resonant active-sterile neutrino conversion 
driven by a pre-existing large lepton asymmetry L = 10~^ — 10^^ which is destroyed dur- 
ing this process. The resonant oscillation between active and sterile neutrinos occurs at a 
temperature 

where £j = 2L^.+J2j=ii ^Vj- The resonant transformation is adiabatic when the active-sterile 
neutrino vacuum mixing is not too small, that is, 

sin^ 2es^ > 10-^ (3) 

for the lepton asymmetry L = 10"'^ — 10~^. Since the resonance temperature of a low energy 
neutrino is higher than that of a high energy neutrino, low energy neutrinos are produced at 
first consuming most of the lepton asymmetry. This makes high energy neutrino conversion 
non-adiabatic and no significant conversion occurs for high energy neutrinos. The produced 
sterile neutrinos are thereby cooler than the active neutrinos and the spectrum is non- 
thermal with {E)/T ^0.7. The free-streaming length of these nonthermal sterile neutrinos 
is 

and the contribution to the matter density today is given by 



In this mechanism, the sterile neutrino production can be confidently calculated only at the 
temperatures below the quark-hadron phase transition temperature, about 150 MeV. This 
translates to the condition 

m.. <23(^^j keV. (6) 

From Eqs. ®-(||), one concludes that about 100 eV to 10 keV sterile neutrinos produced 
through the active-sterile neutrino conversion driven by a lepton asymmetry L ^ 10"^ — 10~^ 
can be a CIDM candidate. 

Our framework is the minimally extended supersymmetric standard model (MSSM) al- 
lowing R-parity and lepton number violating interactions which can explain the neutrino 
masses and mixing implied by the atmospheric and solar neutrino data. In order to es- 
timate the axino mass and mixing with active neutrinos, we need to know the sizes of 
R-parity violating parameters determined from the neutrino data [§,0. The R-parity and 
lepton number violating terms in the MSSM superpotential are 

W = e,^xUH2 + \,kUL,El + K.kLiQjDl , (7) 

where /i is the Higgs mass parameter in the R-parity conserving superpotential, W = fiHiH2. 
As is well known |^ , (active) neutrinos get masses at tree level due to the sneutrino vacuum 
expectation values (VEVs), (L^), misaligned with the bilinear terms e^, as well as at 1-loop 
level through squark and slepton/Higgs exchanges. The tree-level neutrino mass matrix 
takes the form |T7| ; 



^^(M|/Mi/2)c^ where ^ = Jj2^i- Now, the atmospheric neutrino data from Super- 



-^r-e.e.^4' (8) 

where Mz is the Z boson mass, M1/2 is the typical gaugino (photino or zino) mass, and 
C/3 = cos/3. Note that tp = tan/3 = {H2)/{Hi). Here ^i = {L°)/{Hi) - e^ mea- 
sure the misalignment between the sneutrino VEVs and e^. Only one neutrino obtains a 
nonzero mass from the tree contribution (^ and the corresponding eigenvalue is given by 

Kamiokande 0] imply that m^^ ^ y Arriatm ~ 0.05 eV and ^2 ~ ^3 resulting in large mixing 
between the muon and tau neutrino. Therefore, one finds C, ~ ^2,3 and 

ecg ^ 10"^ -^^ -^- . (9) 

^ ^ \Mz J V0.05eVy ^ ^ 

It is important for the later use to recall that nonzero values of C,i arise through renormaliza- 
tion group evolution from the mediation scale Mm to the weak scale and thus can be much 
smaller than the original parameters e^ or A, A' in Eq. (|^). In one step approximation for 



integrating the renormalization group equation (See, for instances, Ref. |jT9|.), one obtains 

(10) 




where hb, Af, are the bottom quark Yukawa couphng and the corresponding soft-parameter 
respectively, and rrij is the slepton mass. Taking fiAi, = rn-andMrn = lO^m^-, one gets 



-2 / ^iC/3 



The 1-loop contributions give nonzero masses to the other two neutrinos, which are much 



smaller than the tree contribution, in particular, in the context of GMSB models [0. Among 
various 1-loop contributions to the second largest neutrino mass m^^, the diagrams with 
interchange of slepton and charged Higgs become important for large tp and thus relevant 



for generating the solar neutrino mass scale. According to the estimation in Ref. |T^, the 
mass ratio m^^/m^^ is given by 



-2 / -^233 




lO-^Tf^ ^ (12) 



under the assumption of the usual hierarchy among the trilinear couplings, that is, those 
involving third generations are larger than the others. The solar neutrino data explained by 
the matter resonant conversion, or the vacuum oscillation [Q requires m^^ ^ JAm'^^i ~ 10~^, 
or 10~^ eV, respectively. Therefore, the matter conversion or vacuum oscillation solution 
can be obtained for tp ~ 50 or 10, respectively, assuming A233 ~ A333. Note also that one 
needs A233, A233333 ~ (10^"^ — 10~^)(mi^2/10^^eV)^/^ for the right solar neutrino mass scale 

iJH- 



Let us now turn to the question why the axino can be CIDM in the context of the R-parity 



violating MSSM with gauge mediated supersymmetry breaking. As discussed in Ref. []T3 
if a superfield S carries a PQ charge xs, there arises the following effective Kahler potential 
for the superfield S and the axion superfield $, which nonlinearly realizes the spontaneously 
broken PQ symmetry below its breaking scale fa, 

i^ = ^($t + $)5t^. (13) 

J a 

Upon supersymmetry breaking, the interaction in Eq. ([13D gives rise to a mass mixing 
between the axino and the fermionic component of S given by 

Fs 
mo ~ Xs— (14) 

J a 

where Fs is the F-term of the field S. Without fine-tuning of some parameters, the value of 
Fs is expected to be of the order of Mg when the field S has the mass M5. As a consequence, 
the axino gets a see-saw suppressed mass, rria ~ rrij^/Ms ~ x|M|//^, that is. 

This shows that the axino has a mass in the right range for being a CIDM sterile neutrino, 
given that fa = (10^° — 10^^) GeV and Ms = (10^ — 10'') GeV when S* is a field in a messenger 
(or a hidden) sector of minimal gauge mediation models [^ . 



A sizable mixing of the axino with active neutrinos arises due to the bihnear terms 
in the superpotential (J^) when the lepton doublets are charged under the PQ symmetry. 
Similarly to Eq. (0), the mixing mass between the axino and the active neutrino Vi is 
given by m^j ^ xlF^/ fa where xl is the PQ charge of lepton doublet L and F/, = eifi{H2). 
Therefore, using the estimation of ej in Eq. (|Tl]), one finds 



keV ^ ^ VSOOGeV/ \ fa J VIO" 

Given nia, the adiabaticity condition 6si ^ nisi/ma ^ 10~^ (|]) puts a lower bound on nisi, 
and an upper bound comes from the fact that the see-saw reduced mass of active neutrinos 
rri^Jma should be smaller than the value m^^ ~ 0.05 eV used in Eq. (j^). Taking roughly 
rn^Jma < 10~^ eV, we get the appropriate range of mixing mass 

10-2 (^) eV < m,, < 3 (^) '" eV . (17) 



VkeV; ~ '*~ VkeVy 

From Eqs. ([I5D and ([T6|), one can find that the sterile and active neutrino mixing mass can 
be obtained in a reasonable range of the parameter space under consideration. 

We have now to comment on the lifetime of our sterile neutrino, the axino, which should 
be stable as a dark matter. As we are working with gauge mediated supersymmetry breaking 
which has low supersymmetry breaking scale, the axino may be able to decay into a gravitino 
and an axion. For the supersymmetry breaking scale of order yF ^ 10^ GeV the gravitino 
mass is ^713/2 ^ 1 keV and thus the axino decay can be allowed kinematically. In this case, 
the lifetime of the axino is found to be ||21[ 



ikevy/ y/F y 

rria I \ 105 GeV J 



10-secP£^ TT^ITTT • (18) 



That is, the axino is stable in the cosmological time scale for all the parameter ranges under 
consideration. 

A nontrivial requirement for the CIDM scenario is a large lepton asymmetry, L = 
10""^ — 10^^, around the nucleosynthesis period. This number is hierarchically larger than 
the baryon asymmetry B ~ lO^^'' required for producing the observed light element abun- 
dances through nucleosynthesis. Therefore, in the absence of a cancellation between L^,^, the 
large lepton asymmetry must be generated at the temperature below the electroweak phase 
transition temperature to avoid the transfer between the lepton asymmetry and the baryon 
asymmetry due to the sphaleron effects. In our framework, since neutrino masses come 
from L violating couplings, they could also be the source of lepton asymmetry. Then the 
most promising way of generating such a lepton asymmetry would be the Dimopoulos-Hall 
mechanism |2^ . 



When a late-time entropy production like thermal inflation is introduced for various 
cosmological reasons [^], the lepton asymmetry can be generated through the late-time 



decay of a (thermal) inflaton followed by lepton number violating superparticle decays 
The lepton asymmetry after the reheat of a late-time thermal inflation is then given by 



Lc^bsL— (19) 

where m^ is the mass of a thermal inflaton 0, Tpt is the reheat temperature after the decay 
of 0, and El is the amount of lepton asymmetry per decay. As a possibihty to get a 
large El, let us consider the case where the stau f is the ordinary LSP which occurs in a 
wide range of the GMSB parameter space [0. As a stau is the LSP among the ordinary 
sparticles, it may decay into a tau and a gravitino (f -^ tG), or into a neutrino and a 
charged lepton (f —>■ vl) through A couplings in Eq. (|^. As we want to have a large lepton 
asymmetry, the latter decay rate has to be larger than the former. In other words, we require 
r(f — > vil) K. m^ |Ap/167r > r(f — *■ tG) ~ m|/167rF^, which puts a bound on A depending 
on the supersymmetry breaking scale yF as follows: 

A > ^^ 1 n-6 ( "^^ ^ ' / lOn^eV \ ' 



In view of generating the solar neutrino mass scale m^^ as given in Eq. (|T^ and the following 
discussions, the above inequality can be fulfilled for reasonable values of a/F. The argument 
goes parallel with the stau decay into two quarks through A' couplings which we neglect 
without loss of generality. As far as Eq. ( |20D is satisfied, the stau mainly decays through R- 
parity and lepton number violating interactions. Out-of-equilibrium condition can be met if 
the reheat temperature is smaller than mf/20 which is the typical supersymmetric particle 
decoupling temperature [||. CP violation can come from complex phases of the leptonic 
charged current in case of Majorana neutrino masses, or from the complex L violating 
trilinear couplings A, A'. Let us now assume for simplicity that the right-handed stau f/j is 
the LSP which has decay modes tr — > z/jCj. Via the tree and one-loop interference depicted 
in FIG. 1, it gives 

where Uij is the CKM type mixing element for the charged leptonic current. The factor 2 
comes from the fact that each produces at least two staus. When for instance A233 is the 
largest coupling, Eq. ( pT]) becomes 

El ~ 2«2lm(f/*2f/33) • (22) 

Note that f/22^33 is complex in general for Majorana neutrinos and its phase unconstrained 
by neutrino oscillation experiments can be of order one. Then the large mixing for at- 
mospheric neutrino oscillations ^ implies Im([/22^33) ~ 1/2- Taking the upper limit of 
Tr Ri rrir /20 which is required for a stau to be out of equilibrium, one finds 

^-?-- (23) 

Therefore, the maximal value of lepton asymmetry in this mechanism is L ~ 10"^. For 
L ~ 10~^, Eq. (^ implies rria ^ 2 keV. From Eq. (|]), we obtain Qa ~ 0.2 for rria ~ 2 keV 
and L = IQ-^ 



An even larger lepton asymmetry can arise through the Affleck-Dine mechanism ||2^ 



with a low reheat temperature. There is an interesting class of models where the lepton 
asymmetry L of order 1 can be produced through the Affleck-Dine mechanism. These are 
supergravity models that possess a so-called Heisenberg symmetry, which include no-scale 
type supergravity models and many string based supergravity models. In these models, 
supersymmetry breaking by the inflationary vacuum energy does not lift flat directions 
at tree level, and one-loop corrections gives small negative mass square for fiat directions 
not involving large Yukawa couplings. After inflation these flat directions generate a large 
baryon and/or lepton asymmetry, typically of order 1, through the Affleck-Dine mechanism 
p5| . For example, we can use the well-known LH2 flat direction for this purpose, with the 
lepton number violation and CP violation arising from nonrenormalizable interactions. In 
this mechanism, it is possible to obtain L ^ 10~^ and hence Qa ~ 1- One difflculty is that 
the large lepton asymmetry may accompany with the large baryon asymmetry. Several ways 
to reduce such large baryon asymmetry to the observed level were discussed in Ref . [^ . As 
we require the reheat temperature Tr to be below the electroweak scale, the decay produced, 
or the regenerated axino population can be completely negligible as can be seen from Eq. (^. 

In conclusion, we have presented a well-motivated candidate for nonthermal sterile neu- 
trino dark matter. The PQ mechanism realized in the supersymmetric standard model to 
resolve the strong CP problem predicts the presence of the axino, the fermionic superpartner 
of the axion. The lightness of the axino can arise when supersymmetry breaking is mediated 
by gauge interactions. We have pointed out that the axino mass can fall naturally into 
the demanded range (0.1 — 10) keV, given the range of the supersymmetry breaking scale, 
Ms ^ (10^ - 10^) GeV, and the PQ symmetry breaking scale, fa ^ (10^° - 10^^) GeV. It has 
also been shown that a sizable mixing between active neutrinos and the axino, 6 > 10~^, 
required for it to be sterile neutrino dark matter, can arise from the (bilinear) R-parity vio- 
lation which is introduced to explain the atmospheric and solar neutrino masses and mixing. 
An 'ad hoc' feature of nonthermal sterile neutrino dark matter scenario is that it needs a 
large lepton asymmetry L ^ 10"*^ — 10"^. We have argued that the lepton asymmetry up 
to L ~ 10~^ can be obtained by a late-time (thermal) inflaton decay followed by R-parity 
violating decays of the ordinary LSP r, given the maximal CP violating phases in the Ma- 
jorana neutrino masses. An even larger lepton asymmetry can result from the Affleck-Dine 
mechanism. The reheat temperature is then required to be below the electroweak scale or a 
few GeV. The scenario under consideration needs better understanding why there is a big 
hierarchy between baryon and lepton asymmetries. 



REFERENCES 



[1] X. Shi and G. M. Fuller, Pliys. Rev. Lett. 82, 2832 (1999), |^ito5-ph/9810076 



[2] E. Gawise and J. Silk, Science 280, 1405 (1998); M. A. K. Gross et al, MNRAS 301, 
81 (1998). 

[3] The Super-Kamiokande Collaboration, Y. Fukuda, et.ai, Phys. Rev. Lett. 81, 1562 
(1998). 

[4] For a recent analysis, see J. N. Bahcall, P. I. Krastev and A. Yu. Smirnov, Phys. Rev. 
D58, 096016 (1998). 

[5] L. Hall and M. Suzuki, Nucl. Phys. B231, 149 (1984); A. S. Joshipura and M 
Nowakowski, Phys. Rev. D51, 2421 (1995); F. M. Borzumati, Y. Grossman, E. Nardi 
and Y. Nir, Phys. Lett. B384, 123 (1996); R. Hempfling, Nucl. Phys. B478, 3 (1996) 
B. de Carlos and P. L. White, Phys. Rev. D54, 3427 (1996); A. Yu. Smirnov and F 
Vissani, Nucl. Phys. B460, 37 (1996); M. Nowakowski and A. Pilaftsis, Nucl. Phys 
B461, 19 (1996); H. P. Nilles and N. Polonsky, Nucl. Phys. B484, 33 (1997); E. Nardi 
Phys. Rev. D55, 5772 (1997). 

[6] M. Drees, et al, Phys. Rev. D57, 5335 (1998); R. Adhikari and G. Omanovic, Phys. 
Rev. D59, 073003 (1999); E. J. Chun, S. K. Kang, C. W. Kim and U. W. Lee, Nucl. 
Phys. B544, 89 (1999); V. Bednyakov, A. Faessler and S. Kovalenko, Phys. Lett. B442, 
203 (1998); A. S. Joshipura and S. K. Vempati, |hep-ph/9808232| ; Otto C.W. Kong, Mod. 



Phys. Lett. A14, 903 (1999); S.Rakshit, G.Bhattacharyya and A. Raychaudhuri, Phys. 
Rev. D59, 091701 (1999); D. E. Kaplan and Ann E. Nelson, |hep-ph/990125^ ; Y. Nir 



and Y. Shadmi, JHEP 9905:023 (1999); A. Datta, B. Mukhopadhyaya and S. Roy, 
|hep-ph/9905549| . 
[7] J. E. Kim, Phys. Rep. 150, 1 (1987). 

[8] E. Kolb and M. S. Turner, The Early Universe (Addison- Wesley, 1990). 
[9] E. J. Chun, A. S. Joshipura and A. Y. Smirnov, Phys. Lett. B357, 608 (1995); Phys. 
Rev. D54, 4654 (1996). 
[10] H. P. Nilles, Phys. Rep. 110, 1 (1984). 

[11] T. Goto and M. Yamaguchi, Phys. Lett. B276, 103 (1992); E. J. Chun, J. E. Kim and 
H. P. Nilles, Phys. Lett. B287, 123 (1992); E. J. Chun and A. Lukas, Phys. Lett. B357, 
43 (1995). 
[12] For a review, see G. F. Giudice and R. Rattazzi, [hep-ph/9801271 . 



[13] E. J. Chun, |hep-ph/ 990 12201 (to appear in Phys. Lett. B). 
[14] K. Rajagopal, M. S. Turner and F. Wilczek, Nucl. Phys. B358, 447 (1991). 
[15] S. Chang and H. B. Kim, Phys. Rev. Lett. 77, 591 (1996). 
[16] L. Covi, J. E. Kim and L. Roszkowski, Phys. Rev. Lett. 82, 4180 (1999). 
[17] M. Nowakowski and A. Pilaftsis, in Ref. §; E. J. Chun and J. S. Lee, |hep-ph/ 98 11201 
(to appear in Phys. Rev. D). 



[18] K. Choi, K. Hwang and E. J. Chun, [hep-ph/9811363| (to appear in Phys. Rev. D); D. 



E. Kaplan and Ann E. Nelson, in Ref. [§. 
[19] E. Nardi, in Ref. [^; A. S. Joshipura and S. K. Vempati, in Ref. |^. 
[20] M. Dine, A. Nelson, Y. Nir and Y. Shirman, Phys. Rev. D53, 1658 (1996). 
[21] E. J. Chun, H. B. Kim and J. E. Kim, Phys. Rev. Lett. 72, 1956 (1994); H. B. Kim and 

J. E. Kim, Nucl. Phys. B433, 421 (1995). 



[22] S. Dimopoulos and L. Hall, Pliys. Lett. B196, 135 (1987). 

[23] D. H. Lyth and E. D. Stewart, Pliys. Rev. D53, 1784 (1996); E. D. Stewart, M. 

Kawasakia and T. Yanagida, Phys. Rev. D54, 6032 (1996); K. Choi, E. J. Chun and 

J. E. Kim, Phys. Lett. B403, 209 (1997); E. J. Chun, D. ComeUi and D. H. Lyth, 

|hep-ph/990328q . 
[24] I. Affleck and M. Dine, Nucl. Phys. B249, 361 (1985). 

[25] M. K. Gaillard, H. Murayama and K. A. Olive, Phys. Lett. B355, 71 (1995). 
[26] B. A. Campbell, M. K. Gaillard, H. Murayama and K. A. Olive, Nucl. Phys. B538, 351 

(1999). 



10 



FIGURES 
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FIG. 1. The diagram for tree and one- loop interference 
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